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Abstract 

Bound and scattering states of a non-central extension of the three- 
dimensional Kepler-Coulomb Hamiltonian are worked out analytically 
within the framework of the potential groups of the problem, 5*0 (7) for 
bound states and 5*0(6, 1) for scattering states. In the latter case, the S 
matrix is calculated by the method of intertwining operators. 



1 Introduction 



> 

in 

(N 

In classical mechanics, the reduced Kepler problem has been known for more 
iy~s . than two centuries [I] to admit seven integrals of motion. These are the total 

^^ ' angular momentum, the Laplace-Runge-Lenz (or Hermann-Bernoulli-Laplace ) 

^^ . vector and the total energy. Since there are two relationships between them 

( see, for example, [5] ) only five of integrals of motion are independent. In 
general, a closed system with N degrees of freedom can have at most 2N — 1 
independent integrals of motion [3]. According to the Liouville theorem, the 
system is completely integrable if it allows N integrals of motion (including the 
C^ ' Hamiltonian ) that are independent and in involution ( i.e. Poisson brackets 

of any two integrals are zero). The system is called superintegrable if there 
exist q, 1 < q < N — 1, additional independent integrals of motion. The cases 
g = 1 and q = N ~ 1 correspond to minimal and maximal superintegrability, 
respectively. In quantum mechanics the definitions of complete integrability and 
superintegrability are same, but Poisson brackets are replaced by commutators. 
The first systematic search for quantum integrable one-particle systems with 
scalar potentials was begun by Smorodinsky and co-workers in (31 [SJ |B] and 
continued by Evans in j7j. It was restricted to the cases when integrals of 
motion are first- or second -order polynomials in the momenta. They found 
all superintegrable systems in two and three dimensions with at most second 
order integrals of motion. It turns out that they possess properties making 



them of special interest: for instance, all these potentials admit the separation 
of variables in several coordinate systems and possess dynamical symmetries 
responsible for the separability of the Schrodinger equation. The history of this 
problem and some results may be found in [S]. 

It is well known that the first quantum study of the hydrogen atom was 
based upon the algebra generated by integrals of motion, before the Schrodinger 
equation was published. Later on, Fock [10] and Bargmann [11] recognized that 
the angular momentum and the Laplace-Runge-Lenz vector generate the Lie 
algebra of SO{A) in the subspace of negative energies and the Lie algebra of 
50(3, 1) in the subspace of positive energies. It was realized that the 'acciden- 
tal' degeneracies, i.e. degeneracies not connected with geometrical SO{3) sym- 
metries of the Hamiltonian, are due to the invariance group 5*0(4). Moreover, 
the separation of variables in parabolic coordinates was related to Laplace- 
Runge-Lenz vector [IT]. Later on, Zwanziger [T^j showed that the algebra of 
50(3, 1) may be used to calculate the Coulomb phase-shifts. Ever since, in- 
variance algebras have been determined for many quantum mechanical systems. 
The best known of these systems are the oscillator [TJ] and the MICZ-Kepler 
svstem [T4l 115) . This is a situation in which the Hamiltonian H of the system 
belongs to the centre of the enveloping algebra of some group G, i.e. 

H - /(O), (1) 

where C is the Casimir operator of the invariance group G. For example, in 
the Coulomb bound-state problem, H = —7^/2(0 -I- 1), where O is a Casimir 
operator of 50(4). 

But it could happen that the Hamiltonian H^, can be related to the Casimir 
operator G as 

^.-/(C)U„, (2) 

where i^^ a subspace occurring in the subgroup reduction and L denotes the 
restriction to i]^. In this case the group G describes the same energy states 
of a family of Hamiltonians H^ with different potential strength. (This is why 
the present group G designated potential group [IH]-) Such an approach was 
proposed by Ghirardi [T7] , who worked it out in detail for the Scarf potential 
[18j . It is similar to the approach of Olshanetsky and Perelomov [TS] [5D] , where 
quantum integrable systems are related to radial part of the Laplace operator 
on homogeneous spaces (i.e. to radial part of Casimir operator of second order) 
of Lie groups. 

Ref. pi] proposed a method that permits purely algebraic calculations of 
5-matrices for the systems whose Hamiltonians are related to the Casimir opera- 
tors G of some Lie group G as ([1]) or ([2]) . Namely, the 5-matrices for the systems 
under consideration are associated with intertwining operators A between Weyl 

equivalent representations U^ and U^ of G as 

S^A (3) 



or 



^k. (4) 



respectively. ( The representations U and U^ have the same Casimir eigenvalues. 
Such representations are called Weyl equivalent.) At this stage we note that the 

operator A is said to intertwine the representations U^ and U^ of the group G 
if the relation 

^C/x (g) = C/x {g)A for all geG (5) 

or 

Adm (6) ^dU^ ib)A for all 6 G g (6) 

holds, where dU^ and dU^ are the corresponding representations of the algebra 
Q of G. Equations ([5]) and © have much restriction power, determining the 
intertwining operator up to a constant. 

The potential group approach has been proven to be useful in variety prob- 
lems in one dimension. Recently, it has been used to describe some potentials 
[22j [23l EH [23 classified in [6] . In Ref. [23] it has been shown that the superpo- 
sition of the Coulomb potential with one barrier term [6] could be related to the 
potential group 50(5). Scattering amplitudes for such system are worked out 
in detail in Ref. [24] by using an intertwining operator [21] between two Weyl- 
equivalent unitary irreducible representations of the SO{5, 1) potential group. 

Subject of the present work will be the simultaneous description of bound 
and scattering states of a quantum mechanical system with Hamiltonian 



1x7^-1+ '' - ^1^ + ^2 - 1/4 sl - 1/4 
2 r 2x^ 2y2 ^ 2z 

written in units h ~ in = \, where Si = 0, 1, 2, . . .. We show that 



^ - -o V - T + o^^ + 0..2 + 9.2 i") 



^ = ^ 9^71 (8) 

2(C+f)^ ^^ 

where C is a Casimir operator of 50(7) (for bound states) or 50(6, 1) (for 
scattering states). 

This system was proved to be minimally superintegrable [6], since four in- 
tegrals of motions were explicitly derived, as a consequence of the separability 
of the related Schrodinger equation in two coordinate systems. But in Ref. [26] 
it has been shown that the classical counterpart of Hamiltonian ([7]) is maxi- 
mally superintegrable, i.e. it admits five independent integrals of motion, in- 
cluding the Hamiltonian: four of them derive from separability of the related 
Hamilton- Jacobi equation in different coordinate systems, but the fifth integral, 
first discussed in Ref. [26], is not connected with separability. Moreover, this 
last integral is quartic in the momenta, while the other three are quadratic, and 
has been rederived in Ref. ;23 ^^ ^^ example of application of a more general 
technique. 

2 General formalism 

Let us start the discussion with the fact that the generators of UIR of SO (7) 
( or SO (6, 1)) are 21 independent Hermitian operators A/^^ — — A/^^ (/i, v = 



1, 2, . . . , 7) which obey the commutation relations 

[Mf,^, Max] = i (5/^o-A^z.A + gyxM^a - g^xM^^ - gi,aM^x) (9) 

where 

g^.u = (+,+,...,+,+) for SO (7) (10) 

5m. = (+,+,...,+,-) for SO (6,1) 

There are three independent Casimir invariants which are identicahy muhi- 
ple of the unit in each UIR. In the case of most degenerate representations, they 
are identically zero, with the exception of the second order Casimir operator 

1 7 
/i,i/— 1 

It is well-known that the most degenerate representation of algebra so(7) 
( so(6, 1) ) can be realized in the Hilbert space spanned by negative-energy 
(positive-energy) states corresponding to fixed eigenvalue of the Coulomb Hamil- 
tonian ijC'°"' in six dimensions 

H^'o^' = ^p'-^, 7>0 (12) 



where a; = {xi,X2, ■ ■ ■ ,xe) e R^, Pj = -«gf-, (j = l,.--,6), x'^ = LLi^*^*' 

p^ = X]i=i PiPi- (We are using units with M = h — 1.) However, in order to be 
able to write the relation ([2|) we introduce the following realization 

Mij = A (x) o (x^pj ~ XjPt) o A"-^ (x) (13) 



2 X'' 



XiP^ ~Pt [x -p) +i-Pi 



^X'^ix), {i,j = l,...,6) 
(14) 



M„ = -MT,^\2h\—^X{x)o 

where 

X{x)=[{xl+xl){xl+xl){xl + xl)]'^' (15) 

and 

h = X{x)o(^ip^-^yx-'{x). (16) 

The generators (|13m4p act in the eigenspace of h equipped with the scalar 
product 

i(l)i,(b2) = J riix)(b2{x) dfi{x), xeR"^ (17) 

where d/x (cc) = A^^ (x) dxidx2 ■ ■ ■ dxQ. 



This representation, of course, is unitarily equivalent to the representation 
constructed in the eigenspace of the Coulomb Hamiltonian H^""^^ in six dimen- 
sions. The unitary mapping W which realizes the equivalence is given by 

^r . i^coui ^ $ = A (a;) *^°"' (18) 

The operators (|13m4|) provide most degenerate representations of SO (7) if 
h is negative definite and of 5*0 (6, 1) if h is positive definite. More precisely, 
they define the most degenerate (symmetric) UIR of SO (7) specified by the 
integer number j = 0, 1, . . .(when h is negative definite) and the most degenerate 
principal series representations of SO (6,1) labelled by the complex number 
j = — I + «p, /9 > (when h is positive definite). If we compute the second- 
order Casimir operator pTj) . it becomes 

Let us consider the reduction corresponding to the group chain G D SO (6) D 
SO (4) X SO (2) D SO (2) X SO (2) x SO (2) , where G is SO (6, 1) or SO (7). 
Then, the basis functions can be characterized by the Casimir operators of the 
chain of groups 

C\r^M) = j-(j+5)|j;/M) (20) 

C^O(6)|_^..^^^ ^ /(? + 4)|j;?M) 

C^°W|j;ZM) = m{m + 2)\j;lM) 
C'^^'^^\j;lM} = sl\j;lM} 
C'"^'^'^^\j;lM) = sl\j;lM) 

where Af is a collective index (m, si, S2, S3) and 

C.so(6) ^ i ^ m2, C^oW = - E M^j^C^""^^^' = Af22,C^o(2)2 = M|4,C^o(2)3 ^ ^_ 

(21) 
According to this, we introduce in place of xi, 2:2, . . . , xg the variables r, 0, 99, ai, a2, 0^3 
via Xi = rui with 

(22) 



rii 


= smt/suKpsmai, 


?^2 ^ 


= smf^suKpcosai 


ri3 


= sin6'cos(^sina2, 


^4 


= sin 9 cos (/3 cos a2 


n-5 


— cos0sina3. 


n-6 


= cos 6 cos a3 



where < r < 00, < 6*, (/? < ^ and < ai, 02, a3 < 27r. If we compute the 
operator 7^/ (C -I- ^) for this parametrization, it becomes 



2 

7 


92 2 9 1 

Qj.2 J. Qj, j,2 

1 
1 


\siii6de"'^ 86 ' sin2 9(^2J 

a 92 \ 1 /I 92 ^ 

III 1 


, 1 


C9 


C+f 


/I 92 
U + 9a2 




r^ sin^ sin^ ^p 


^4 ' dal) ' r2 sin2 6* cos2 (^ V4 ' ^a^y 


' ' r2 cos2 ' 



Let 'HsiS2S3 be a subspace spanned by |j; IM) with fixed si, S2 and S3. Thus, 
the operator (P5|) restricted to this subspace becomes a differential operator in 
r, and Lp\ it turns out that 



7' 



C 



25 
4 



, 1/4- s? , 1/4- si , 1/4- si 



r^ sin sin if r^ sin cos^ ip 
with Si = 0, ±1, ±2 . . ., where we have used that 

daf 
Hence, the Hamiltonian 



H = -V-^+ ;'^';^/', (25) 

2 r 2r2 sin^ 6* sin V 

, ^i-1/4 , ^i-1/4 



2r2 sin^ 9 cos^ (^ 2r2 cos^ 9 
can be described in terms of the potential groups SO (7) and SO (6, 1) since 



7' 



C+f 

as mentioned in Section 1 (formula (|H1))- (Due to the symmetry st — >■ — s^ in 
the Hamiltonian (j25l) , without loss of generality, we may assume that si, S2 and 
S3 are non- negative integers.) Note that the 5*0(2) subgroups are related to 
potential strength. 

At this point, it is worthwhile pointing out that Hamiltonian (j25p does not 
contain the pure Coulomb potential as a particular case, within the framework 
of the 50(7) and 5*0(6,1) symmetries considered in the present work. In 
order to restore it, it is necessary to resort to larger symmetry groups, for 
example, SO (10) and SO (9, 1) and use the decomposition chain G D SO (9) D 
5*0 (6) X SO (3) D 5*0 (3) x SO (3) x 5*0 (3) , where now the 5*0(3) subgroups 
are related to potential strength. 

Here again, use is made of polar coordinates 

X = (sin0siniy9ei, sin cos (/9e2 , coshes) 

where x £ R^ , Ci = (sin ai sin Pi , sin ai cos (5i , cos ai),i — 1,2,3. Then, a proce- 
dure similar to that described above would lead to the Hamiltonian 

^^_1^2_7, ^l(^l + l) , ^2(^2 + 1) ,^3(^3 + 1) 



2 r 2r2 sin^ 9 sir? ip 2r^ sin^ 9 cos^ ip 2r'^ cos^ 



where li {i ~ 1,2,3) are integer and are allowed to take the null value, thus 
restoring the pure Coulomb potential. 
Finally, we note that the operators 

«2 _ 1 «2 _ 1 «2 _ 1 

r _ T 2 I '^1 4 I ^2 4 , f3 4 /qra 

^1 - ■'^ + ..:, 2/1..:. 2,, + ,.:„2/i_„9., + 2 ^ V^^) 



h = L 



s\v? 6 s\v? ip sin2 cos2 (/j 

s2 _ 1 o2 _ 1 

2 I *1 4 I *2 4 



sin2 (^ cos2 (^ 



where L2 and i2 ^j-g the square of "angular momentum" and of its projection 
on the third axis, commute with the Hamiltonian. These integrals of motion 
are related to the Casimir operators of SO (6) and its SO (4) subgroup in the 
sense that 



where 

(jSOCo) 



/i = c^o{&^ 


, J2 = C^O(4) 


■HsiS2S3 


(27) 


( I d d 1 92 \ 






1 /I 92 \ 1 

1 


/I 92 \ 

1 


1 /I 92 

1 


sin^ 9 SIT? <^ \^ da\ J sin^ 9 cos^ 


if, V4 ' dalJ 


COS^^ \A ' 9^3 


4)_ 52 

9(/32 s, 


1 /I ^ 92 \ 


1 /I , 92 \ 




n2(^V4 t?"iy cos2(^y4 da^J 





and 



2.1 Bound states 

The bound state spectrum can now be easily obtained if we note that the eigen- 
value of the Casimir operator C of the potential groups SO (7) is j {j + 5). We 
then find 

E = ^ ^ (28) 

20 + 1)' 

where j = si+ S2 + S3 + 2fci + 2fc2 + n, (fci, ^2, n = 0, 1, 2, . . .). It is easy to 
check that states (|28|) have degeneracy 2 ^ where d — [ -'^''^^''^^''^ j + 1, and 
[q\ is the largest integer less than or equal to q. 

We give for reference the expression of the bound-state wave functions 

^(x)=7^,^(r)3^^M(0,^), (29) 

where TZji (r) is the radial part of the wave function, while yiM (^7 95) is the 
angular part of it : 

7^,, (r) = c«'+ie-tL2^+4 (^) ^ ^ ^ 2^r / (j + ^) (30) 



with n = j — I (n = 0, 1, 2, . . .), 



27 

J + l 



r{j-i + i) 



2{j + l)rij + i + 5) 



and 



yij^,j (e, ip)^x sin™+^ e cos"3+^ e sin"i+^ (^ cos^^^+s (^ 

j{m.+l,Si) (^^^on\ pC^'i^'Sa) 
fci 



^p,„.-r.,«; (cos26l)Pf ^^""^ (C0S2(^) 



(31) 



(32) 



with 2ki ^ I — m — S3, 2k2 — m -- si — S2 {ki, fc2 = 0, 1, 2, , . . .) and 

'r{^{l + m + S3 + 4)) r (^ (/ - m - S3 + 2)) T (^ (m + -si + S2 + 2)) 

r (i (; + TO - S3 + 4)) r (i (/ - to + S3 + 2)) r (i (to + si - S2 + 2)) 

- S2 + 2)) 



r(i(TO 



r(i(TO-Sl+S2+2)) 



{21 + 4) (2to + 2) 



(33) 



5(",/3) 



Here, L" and PH are Laguerre and Jacobi polynomials, respectively. 
It is worth noting that (see Appendix ) the yiM functions are related to 5- 
dimensional spherical harmonics Yim (n) (see Section 10.5 of [28) ) in polyspher- 
ical coordinates, while TZji is related to the radial part of the 6-diinensional 
Coulomb wave function [29] as 

n,i (r) = r^TZfr' (r) 

2.2 Scattering states 

Once the group structure of the problem has been recognized, the associated S 
matrix can be computed by using Eqs. I^S^. This requires knowledge of ma- 
trices {I'M' \A\ IM) that intertwine Weyl-equivalent representations of SO (6, 1) 
in the bases corresponding to the SO (6, 1) D SO (6) D SO (4) x SO (2) D 
SO (2) X SO (2) X SO (2) reduction. One has (see Appendix) 



where 



{1'M'\A\IM)^Ai5iv5mm' 
T{^ + ip + l) 



A, 



According to this, we have 

S {9, ip; 9', ^') = J2 ^'^'^^ (^' ^) y^M {0\ f') 



IM 



Thus, the scattering amplitude, / {9, ip; 9' , ip'), is defined by 

/ {9, cp; 9', ^') = — E (^' - 1) ^"^ (^' ^) yiM (^'' V) 



(34) 
(35) 

(36) 

(37) 



Since 

Y, yiM {0, ^) ym [9', if') = S (cos - COS 0') 5 [ip - ip') 

IM 

we can omit unity in the brackets of formula (j37p when 6 ^ 6' , Lp ^ ip\ leaving 
/ (0, ^; 0', ^') = — E ^'^"^ (^' ^) yiM (0', V') ■ (38) 

^ IM 

Moreover, formulas (|49p . (l5^ and pS)) imply the following integral represen- 
tation of the scattering amplitude 



2-K 27r 27r 

f {9,ip;0' ,(f') — ^rjvb / / (1 — a sin sin 6*' — cos cos 0' cos as) 





5-ip 



where 
and 



X exp {—isiai — is2a2 — ^s^a-^) daida2da^ (39) 

a — sin Lp sin ip' cos ai + cos (^3 cos ip' cos a2 (40) 

& = sin 6 sin sin Lp sin (^9 (41) 



3 Conclusions and outlook 

We have shown in the present work, based on the potential group approach, 
how a non-central extension of the Coulomb Hamiltonian, considered in the 
literature as an example of maximal superintegrability, can be worked out in a 
fully analytic way, with bound states described by most degenerate represen- 
tations of SO (7) and scattering states by most degenerate representations of 
5*0(6, 1). The subfamily of the generalized Coulomb problem described in the 
present work does not include the pure Coulomb potential: in order to restore 
it, the symmetries could be enlarged to 5*0(10) and 50(9, 1), respectively. 

The generation of solvable non-central potentials via the potential group ap- 
proach is quite general and not limited to the orthogonal and pseudo-orthogonal 
groups of interest to the Coulomb problem. An example of a non-central ex- 
tension of the harmonic oscillator with U (4) symmetry has been discussed in 
Ref. [35] , while a non-central extension of the null potential with E (4) symme- 
try has been worked out in Ref. [25, . Other cases of physical interest with more 
complicated symmetries will be considered for future work. 

4 Appendix: Calculation of the matrix elements 
of A 

Here we calculate the matrix elements of A which intertwine Weyl-equivalent 
representations of 50(6, 1) or so (6, 1) in the bases corresponding to 50(6, 1) D 



SO{6) D SO (4) X SO (2) D SO (2) x SO (2) x S'O (2) reduction. We find it 
expedient to use, for tfiis purpose, equation (|6]). 

We shall start with the fact that the most degenerate principal series repre- 
sentations of SO{6, 1) can be realized on £2 (<5'^) (see Section 9.2.1 of [28] ) 

C/,(5)/(n) = Kr7K) , neS' (42) 

where 

6 v^6 -1 , 

E-i , / ^ Z^i=i 9ki "-« + 9k7 

9n n, + 377, {rig)^ = \ ^_\ ■ 

The operator A defined by 



{Af){n)^ J K{n,n')f{n')dn' (43) 

intertwines representations j and — 5 — j, if 

K{n„n'^) = (0.,)^+^ (c.;)'+^if(n,n') . (44) 



The kernel, K , is uniquely determined by Eq. (I44p up to a constant and is 
given by 

K{n, n') =ri{l-n- n')"^"^ • (45) 

with 

7? = 2-^+^"^^^^ (46) 

TraF {—ip) 

With this factor the operator A becomes unitary for j = —^ + ip (see equation 

®)-. . 

Taking into account the fact that 5-dimensional spherical harmonics Yi]\j of 

degree I [28] forms a bases in £2 (S^), corresponding to above reduction, we 

have the following integral representation for the matrix elements of A 



{l'M'\ A \IM) = / /C (n, n') Y*^' {n') Yim {n) dndn' . (47) 

where dn — sin 6 cos 6 sin ip cos Lpd9d<fda3da2dai for n as in (1221) and 

3 1 

YiM (n) = yiM (9, p) n -=e^^^-"^- . (48) 

By using the expansion 

^ (1 - n . n!ri-'^ ^ 1 g (. + 2) l%^''^/\ cl {n ■ n') , (49) 

2^ ^ r (2 - «p + j^) 

we have 

{1'M'\A\IM)=Ai5u'5mm' (50) 



10 



with 

A. = £i|±^ (51) 

In arriving at equation (j50p we have used the addition formula 

^' ("•"') = -^ E ^-^M (n) r;M (^') (52) 

References 

fll Goldstein H 1975 Am. J. Phys. 43 737; Goldstein H 1976 Am. J. Phys. 44 
1123; 

Goldstein H, Poole C and Safko J, Classical Mechanics, Addison- Wesley, 
2001. 



[6: 



[9 

[lo; 
[11 

[12 

[13; 

[14 

[is: 

[16; 

[17; 
[is: 



Landau L and Lifshitz E M, Mechanics, 3rd ed, Pergamon, Oxford, 1976. 

Fris I, Mandrosov V, Smorodinsky YaA, Uhlir M and Winternitz P 1965 
Phys. Lett. 16, 354. 

Winternitz P, Smorodinsky YaA,, Mandrosov V, Uhlir M and Fris I 1967 
Sov. J. Nucl. Phys. 4, 444 

Makarov AA, Smorodinsky YaA, Valiev K and Winternitz P 1967 Nuovo 
Cimento A 52 1061. 

Evans NW 1990 Phys. Rev. A 41 5666. 

Winternitz P 2009 Physics of Atomic Nuclei 72, 875. 

Pauh W 1926 Z. Phys. 36 336. 

Fock V 1935 Z Physik 98 145. 

Bargmann V 1936 Z Physik 99 576. 

Zwanziger D 1967 J. Math. Phys. 8 1858 

Jauch JM and Hill EL 1940 Phys. Rev. 57 641. 

Zwanziger D 1968 Phys. Rev. 176 1480 

Mcintosh HV and Cizveros A 1970 J. Math. Phys. 11 896. 

Alhassid Y, Giirsey F and lachello F 1983 Ann. Phys. 148 346. 

Ghirardi GC 1972 Nuovo Cimento A 10 97. 

Scarf F L 1958 Phys. Rev. 112 1137. 



11 



[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 
[27] 

[28] 

[29] 



Olshanctsky MA and AM Pcrclomov 1977 Lett. Math. Phys. 2 7. 

Olshanctsky MA and AM Perelomov 1983 Phys. Rep. 94 313. 

Kerimov GA 1998 Phys. Rev. Lett. 80 2976. 

Kerimov GA 2006 Phys. Lett. A 358 176. 

Kerimov GA 2006 J. Phys. A: Math. Gen. 39 1183. 

Kerimov GA 2007 J. Phys. A: Math. Theor. 40 7297. 

Kerimov GA 2007 J. Phys. A: Math. Theor. 40 11607. 

Verrier PE and Evans NW 2008 J. Math. Phys. 49 022902. 

Rodriguez MA, Tempesta P and Winternitz P 2009 J. Phys. Conf. Ser. 175 
012013. 

Vilenkin NJa and Khmyk AU, Representation of Lie Groups and Special 
Functions Vol.2, Kluwer Academic Pubhshers, Dordrecht, 1993. 

Nieto M M 1979 Am. J. Phys. 47 1067 



12 



